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Abstract. We present dilepton spectra from p+p, d+p and p+Nb reactions at SIS energies, which were 
simulated with the GiBUU transport model in a resonance model approach. These spectra are compared to 
the data published by the HADES and DLS collaborations. It is shown that the p spectral function includes 
non-trivial effects already in elementary reactions, due to production via baryon resonances, which can 
yield large contributions to the dilepton spectrum. Dilepton spectra from nuclear reactions in the energy 
range of the HADES experiment are thus found to be sensitive also to properties of nucleon resonances in 
the nuclear medium. 

PACS. PACS-key discribing text of that key - PACS-key discribing text of that key 



o 

=5 



> 

in 
in 

en 

en 

o 

(N 



X 



1 Introduction 

While the vacuum properties of most hadrons are known 
to reasonable accuracy nowadays, it is a heavily debated 
question how these properties change inside nuclear mat- 
ter. In particular, various theoretical predictions regarding 
the in-medium properties of the light vector mesons have 
been suggested. For recent reviews on in-medium effects, 
see [Tppj. 

Among these expected in-medium effects, a so-called "col- 
lisional broadening" of the meson spectral function, due to 
collisions with the hadronic medium, is expected. A second 
class of predictions claims that the vector-meson masses 
are shifted in the medium due to the partial restoration of 
chiral symmetry 0]. QCD sum rules can constrain these 
effects, but do not provide definitive predictions [5j. 
The more prominent hadronic decay modes of the vector 
mesons are unfavorable for studying in-medium effects, 
since they are affected by strong final-state interactions 
with the hadronic medium - in contrast to the rare dilep- 
ton decay modes. As the leptons only interact electromag- 
netically, they are ideally suited to carry the in-medium 
information outside to the detector, nearly undisturbed 
by the hadronic medium. 

Dilepton spectra from nuclear reactions with elementary 
projectiles have been studied for example with the CLAS 
detector at JLAB, where photons with energies of a few 
GeV interact with nuclei [6], or by the E325 experiment 
at KEK, where 12 GeV protons were used as projectiles 
[7 . On the side of the hadronic decays, most notably u —> 
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7r°7 is being investigated by the CBELSA/TAPS group in 
photon-induced reactions 0[!5]- 

Early measurements of dilepton spectra from heavy-ion 
collisions in the low-energy regime were conducted by the 
DLS collaboration [10] . showing an excess over the ex- 
pected yield. A similar excess was also observed in ex- 
periments at higher energies [TTlll2| . where it could be 
attributed to an in-medium broadening of spectral func- 
tions [13, 14, 15] . For the DLS data such in-medium effects 
never provided a convincing explanation - a problem that 
was soon known as the "DLS puzzle" [15111711151115] . 

More recently, the HADES collaboration at GSI has set 
up an ambitious program for measuring dilepton spectra 
from p+p, p+A and A+A reactions [2lJ l l2"T P 2 l l2^ 1 124 l l23 ] . 
in order to systematically check the old DLS data with 
improved statistics and to finally resolve the DLS puzzle. 
Up to now this endeavor has fully confirmed the validity of 
the DLS data and shifted the puzzle into the theory sector. 
It is clear that a detailed understanding of the elementary 
reactions is the most important prerequisite for explaining 
the heavy-ion data. 



In this paper, we apply the Giefien Boltzmann-Uehling- 
Uhlcnbeck transport model (GiBUU) 26 to the elemen- 
tary reactions (nucleon-nucleon and proton-nucleus) stud- 
ied by the HADES collaboration. We use GiBUU to gen- 
erate dilepton events and pass them through the HADES 
acceptance filter, in order to compare our calculations di- 
rectly to the experimental data measured by HADES. 
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2 The GiBUU transport model 

Our tool for the numerical simulation of dilepton spectra 
is the GiBUU hadronic transport model, which provides 
a unified framework for various types of elementary re- 
actions on nuclei as well as heavy-ion collisions [26,27 . 
This model takes care of the correct transport-theoretical 
description of the hadronic degrees of freedom in nuclear 
reactions, including the propagation, elastic and inelastic 
collisions and decays of particles. 

In GiBUU the spectral one-particle phase-space distribu- 
tions, F(x,p), of all particles are obtained by solving the 
coupled Kadanoff-Baym equations [28 for each particle 
species in their gradient-expanded form |29j 



VF(x,p)-tr{rfReS™ t (x,p)},=C(x,p) , (1) 



with 



VF = { Po -H,F} 



P b 



(2) 



Here {. . .} p b denotes a Poisson bracket. In the so-called 
backflow term (second term on the left-hand side in (fil), 
that is essential for off-shell transport, f(x,p) is the phase- 
space density related to F by 



F(x,p) = 2wgf(x,p)A(x,p) 



(3) 



where A(x,p) is the spectral function of the particle^] and 
g is the spin-degeneracy factor. The quantity _T in the 
backflow term is the width of the spectral function, and 
S ret (x,p) denotes the retarded Green's function. Off-shell 
transport is thus included and leads to the correct asymp- 
totic spectral functions of particles when they leave the 
nucleus. The expression C(x,p) on the right-hand side 
of (TTT) denotes the collision term that couples all parti- 
cle species; it contains both a gain and a loss term. For 
a short derivation of this transport equation and further 
details we refer the reader to [26 . In order to solve the 
BUU equation numerically, we rely on the test-particle 
ansatz. Here the phase-space densities are approximated 
by a large number of test particles, each represented by a 
(^-distribution in coordinate and momentum space. 
The collision term contains all sorts of scattering and de- 
cay processes: elastic and inelastic two-body collisions, de- 
cays of unstable resonances and even three-body collisions. 
The two-body part of the collision term is separated into 
two different regimes in terms of the available energy, \fs: 
a resonance model description at low energies and the Py- 
thia string model at high energies. 

For baryon-baryon collisions, the transition between the 
two is usually performed at \fs = 2.6 GeV. There is a 
small window around this border (±0.2 GeV), where both 
models are merged linearly into each other in order to 
ensure a smooth transition. For meson-baryon collisions, 
the transition region lies at \fs = 2.2 ± 0.2 GeV. 
Unfortunately, the transition region in this default GiBUU 
prescription lies right inside the range of energies used for 

1 A is normalized as J °° A(x,p)dpo = 1. 



the HADES experiment. However, we think that it is im- 
portant to describe all HADES spectra with one consistent 
model. In this paper we therefore explore the possibility 
of pushing the transition region up to higher energies and 
using an extended resonance model for all reactions mea- 
sured by HADES. 

In the high-energy regime the GiBUU collision term relies 
on the Monte Carlo event generator Pythia (v6.4) [30] 
[3Tj . which is based on the Lund string model. Although 
Pythia clearly has its strengths at higher energies (tens 
to hundreds of GeV), it is used in GiBUU down to ener- 
gies of a few GeV. This works surprisingly well, as has 
recently been demonstrated for example by GiBUU's suc- 
cessful description of pion data measured by the HARP 
collaboration [32] . 

Despite this good description of pion obscrvables in the 
few-GeV energy regime, it turned out that the HADES 
dilepton data for p+p collisions at 3.5 GeV pose a some- 
what greater challenge for Pythia [33J . Most prominently, 
the vector-meson production is strongly overestimated by 
the default Pythia parameters, and also the intrinsic px 
distribution needs to be adjusted slightly to reproduce the 
HADES p T spectra. 

Since a resonance description should in principle be appli- 
cable in the energy regime probed by the HADES exper- 
iment (y/F < 3.5 GeV), we try in the following to set up 
such a description as an alternative to the string model 
approach. 



3 The resonance-model approach 

The low-energy part of the nucleon-nucleon collision term 
is given by a resonance model based on the Teis analy- 
sis [34], in which all collision cross sections are assumed 
to be dominated by the excitation of baryon resonances. 
The GiBUU model currently contains around 30 nucleon 
resonances, for a complete list see [26]. However, only the 
subset used in the Teis analysis is actually being popu- 
lated in NN collisions, see table!!] The properties (masses, 
widths and branching ratios) of all the resonances are 
taken from the partial- wave analysis of Manley [35] . All of 
these states, except for the P33(1600), are not only found 
in the Manley analysis, but have been confirmed, e.g., by 
the more recent analysis of Arndt et al. [36] and received 
a four-star rating from the PDG [37]. We note already 
here that some of the branching ratios which are impor- 
tant for the present study, in particular those for decay 
into pN and ujN are not very well known and still under 
experimental investigation J38tl39j. 

We use all the resonance parameters and branching ratios 
exactly as given by Manley, with one exception: The pA 
decay channels are introduced by Manley only in order to 
account for missing inelasticities, which are not covered 
by one- and two-pion final states. In that sense, Manley 
has no real evidence for the pA final state in particular, 
but just uses this decay channel to account for the left- 
over strength. Therefore we take the freedom to replace 
the pA decays by a A, in order to avoid an overestimation 
of the p-meson production. The influence of Manley's pA 
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Table 1. Resonance parameters according to Manley [35] (columns 2-4), together with matrix elements for production in pp 
collisions (columns 5 and 6) and branching ratios of the resonance decay modes (columns 7-13). Subscripts indicate the relative 
angular momentum of the outgoing particles in the respective decay channel. 



decay channels on dilepton spectra was already discussed 
in [3D] for the case of pion-induced reactions. The dilepton 
spectra actually give a hint that the needed Sir inelasticity 
might not be in the pA, but instead in some other channel, 
as e.g. a A. 

Also the width parametrizations are taken from the Man- 
ley analysis, where the partial widths for, e.g., A — > ttN 
and p — > 7T7T are parametrized according to 



r(m) = r a 



m 
m 



( lo 



3 2 
% 



A 2 



A 2 ' 



(4) 



Here too is the mother particle's pole mass, to is its off- 
shell mass, Jo is the on-shell width (at to = too); q de- 
notes the final-state center-of-mass momentum for mass 
to, while go is the same quantitiy for mass too, and A = 
1/ R = 1 fm _ can be viewed as a cutoff-parameter. It has 
been shown in [51], that eq. Q gives a good description 
of the experimental phase shifts in 7T7t and ttN scattering. 
For the detailed treatment of the other decay channels, we 
refer to chapter 3.3.1 of [26 . 

The resonance model used in this work is based on the Teis 
model, but modifies and extends it in several aspects. We 
take into account the following nucleon-nucleon scattering 
channels: 

1. NN -> NN 

2. NN -> NA, 

3. NN -> NN*, NA*, 

4. NN ^ AA, 

5. NN -> AN*, AA*, 

6. NN -> NNir (non-res. BG) 

7. NN -> NNu, NNttuj, NN<j) (non-res.), 

8. NN -> BYK (with B = N,A;Y = A, £). 

For the elastic cross sections (first item), we rely on the 
parametrizations by Cugnon et al. [42j (for beam momenta 
below p uh w 2.776 GeV) and the PDG [33] (above). For 
details see also l26l. 



The single-resonance excitation channels (items 2 and 3) 
were already included in the Teis analysis. While the NA 
channel is treated by an OBE model according to Dmitriev 
et al. [44] . the higher resonances are produced in a pure 
phase-space approach with constant matrix elements, 



&NN^NR 



Ci \M 



NR\ 



PiS 16ir 



dnA R (p)p F {p)- (5) 



Here, pi and pp denote the center-of-mass momenta in the 
initial and final state, respectively. The matrix elements, 
Mnr, have previously been fitted by Teis to exclusive 
meson production (n, 2ir, r\ and p). Our values are listed 
in tab. [TJ Ar denotes the resonance spectral function, 



A R (p) = - 



^TRip) 



{p 2 -MlY + p 2 rl{p) 



(6) 



In principle all production channels are assumed to be 
isospin-symmetric, with the Clebsch-Gordan factors, Ci, 
resulting from this symmetry. The only exception from 
this isospin symmetry is the 5n(1535) resonance: The ex- 
clusive 7] production, which is assumed to proceed exclu- 
sively via this resonance, is known to be significantly larger 
for pn than for pp [45] , therefore we use 

\M.pn-^.NN*(1535)\ = 6-5 • \M pp ^NN- (1535) I ■ (7) 

Note that while the Sii(1535) is known to dominate the r\ 
production in pp at low energies, there may of course be 
other contributions [4"fj] . 

The singlc-pion production cross section can not be de- 
scribed satisfactorily by resonance contributions alone, and 
one has to add a non-resonant background term [2"6ll3"¥] 
(slightly refitted here), whose largest contributions appear 
on the left-hand shoulder of the NA peak. 
Most of the resonance-production matrix elements are ad- 
opted from Teis. However, me make a few modifications. In 
particular we reduce the contributions of the Di5(1675), 
P3i(1910) and P33(1600), which were extremely large in 
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Fig. 1. (Color online) Inelastic pp and pn cross sections in the resonance model. The data points shown here have been obtained 
by subtracting the parametrized elastic cross section from the total cross section data |37| . 



the Teis analysis, in favor of the Pn(1440) and double- A 
contributions. This gives an improved threshold behavior 
of the 2-7T production channels (in line with the analysis 
of Cao et al. 5Z])j as weu as a better agreement with the 
inelastic nucleon-nucleon cross sections, cf. fig. [I] 
Furthermore, we add another isospin-asymmetry factor for 
the Pu(1440) state: 



\M 



pn^-NN* (1440) I 



= 2-\M 



pp-yNN* (1440)1 



(8) 



This is done in order to improve the agreement with the 

np inelastic cross section data, which would otherwise be 

underestimated significantly. 

Also the double-resonance production (items 4 and 5), 

which in the Teis model was limited to A A, is performed in 

a phase-space approach, analogous to the single-resonance 

excitation: 



PNN^AR 



Ci \M 



AR\ 



PiS 



16tt 



/ d/iid^A^MiMfl^WOi:^)- 



(9) 



Here one integrates over the spectral functions of both 
resonances (/Xi,2 being their masses). 
In the Teis analysis, the production mechanisms were re- 
stricted to NN -> NR and AW -> AA, so the obvious 
extension candidate would be general double-resonance 
excitation channels (AW —> R1R2). The channels taken 



into account by Teis were fitted to single- and double- 
pion production data. Therefore his model is only guar- 
anteed to work in the low-energy region. At higher en- 
ergies, the model starts to fail, since the more inclusive 
multi-meson final states are not included. If we want to 
describe AW collisions in the HADES energy regime of 
yfs « 2.4 — 3.2 GeV with a resonance model, we clearly 
need to extend the Teis approach. We do this by restrict- 
ing ourselves to the same set of resonances (cf. tab. fl]), 
but extending the production mechanisms. 

Since Teis already describes the exclusive n, 2n, p and r\ 
production, what is missing are channels like e.g. nr\, np, 
37T, 2i], 2p, etc. Unfortunately there are almost no expe- 
rimental data available for these channels. We thus have 
to rely on the cross sections obtained from Pythia as an 
estimate to fix these channels. According to Pythia, the 
inclusive p and r\ production is in fact dominated by the 
channels np and nr/, respectively, at the highest HADES 
energy of y/F « 3.2 GeV. Therefore we concentrate on 
these two for now, and neglect all others. Our strategy to 
satisfy these channels relies on double-resonance excita- 
tion, NN —t AR, where the A decays into ttN, while the 
other resonance R will be one with an r/N or pN decay 
channel, so that we end up with a nrj or np final state 
(note that we do not include cascade decays of single res- 
onances into nrjN, as treated for example in |48j . since 
our model misses the corresponding decays modes, such 
as TjA). We add three new classes of production channels: 

i) NN -> Z\Sii(1535) (-> NNnrj), 
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ii) NN -» AN* (-> NNnp), 

N* = D 18 (1520), S u (1650), F 16 (1680), P 13 (1720), 

iii) NN -> Z\Z\* (-> NNnp), 

A* = S 31 (1620),D 33 (1700),F 35 {1905). 

For each of these we need one new parameter, namely the 
matrix elements, |A^i| 2 /16tt = 60mbGeV 2 , |A^ 2 | 2 /16tt = 
12mbGeV 2 and |_M 3 | 2 /167r = 21mbGeV 2 , as listed in 
tab. T] As noted before, we fix the matrix elements to 
roughly fit the Pythia cross sections for ttt] and wp pro- 
duction (with further constraints from the total pp cross 
section as well as the HADES dilepton data). As for the 
exclusive production, we assume that the rj meson is pro- 
duced exclusively via the 511(1535), while the p produc- 
tion proceeds via a number of N* and A* resonances. 
It should be noted that the new channels will not affect 
the exclusive meson production fitted by Teis, except for 
the 27r channel, which gets minor contributions from these 
channels. 

The production of u and <fi mesons is not carried out via 
baryonic resonances in our model (although a coupling of 
the u) to nucleon resonances has been reported for example 
in [49,50,51,38 ). Instead, their production cross sections 
are parametrized in a phenomenological manner |52) : 



a(pp — > ppV) — a(l — x) x c , with x = Sq/s. (10) 

Here sq = (2itin + my) 2 is the threshold energy, and the 
parameters a, b and c are listed in table [2j We use this 
parametrization not only for exclusive w and <fi produc- 
tion, but also for NN — > NNirui. Since there are no data 
available for this channel, we fitted its parameters to the 
Pythia results. 





v^ [GeV] 


a [mb] 


b 


c 


Ref. 


TTU) 



2.658 
2.796 
2.895 


5.3 
1.0 
0.01 


2.3 

1.5 
1.26 


2.4 
1.1 
1.66 





Table 2. Parameters for vector-meson production. 



As seen in fig. [II we achieve a good agreement with data 
for the inelastic pp cross section up to about yfs = 3.5 GeV. 
At higher energies 37r and 4tt production becomes impor- 
tant, which is underestimated by our model (and other 
channels which we miss completely). In the np cross sec- 
tion there are minor deviations, and unfortunately also 
the quality of the data is not quite as good as for pp. 



4 Dilepton decays and form factors 

In the GiBUU model the following dilepton decay modes 
are taken into account: 



• direct decays, asF-> e + e , 
with V = p°, ui,<f> or 77 — > < 



Dalitz decays, as P 
or ui — > 7r°e + e~ c 



+ e 7 with P = 7T° 
A -> Ne + e~ . 



V 



Most of them are treated similarly as in [JT] . The leptonic 
decay widths of the vector mesons are taken under the 
assumption of strict vector-meson dominance (VMD), 



r v 



5 -(m) = c 'v' 



,3 ' 



(11) 



where p, is the meson's off-shell mass, my is the pole mass, 
and the constants Cy are listed in table [3] (taken from 
[37]). Although the physical threshold of the dileptonic 
decay channels of course lies at 2m e , contributions of p 
mesons below m = Im^ are frequently neglected in trans- 
port simulations. The reason for this artificial threshold is 
purely numerical: The p spectral function has a sharp drop 
at the 27T threshold, and it is numerically very difficult to 
populate the spectral function below this threshold, where 
it is almost vanishing. Here we make additional numerical 
efforts to include the contribution of p mesons below the 
27T threshold, since it can give significant contributions to 
the total dilepton spectrum for certain reactions. 



V m ¥ (MeV) r ee (keV) 



P 



775.49 

782.65 

1019.455 



7.04 
0.60 
1.27 



Cy — r^/mv 



9.078- lO^ 15 " 

7.666 ■ 10~ 7 

1.246 ■ 10~ 6 



Table 3. Dilepton-decay constants for V 



While the direct decay of the 77 meson into a p + pT pair 
has been observed, for the corresponding e + e~ decay only 
an upper limit of BR(r7 — > e + e~) < 2.7T0~ 5 is known |55j . 
In fact this limit has been pushed down to 4.9 • 10~ 6 lately 
using HADES dilepton data (21] • However, the theoretical 
expectation from helicity suppression is still four orders 
of magnitude lower [56 . The absence of any 77 peak in 
the measured spectra allows us to conclude that the true 
branching ratio must be significantly lower than the upper 
limit just mentioned |57j . Therefore we do not include the 
rj — > e + e~ decay in our analysis. 

The Dalitz decays of the pseudoscalar mesons, P = tt°,7], 
are treated via the parametrization 58J, 



dp 3-7T p 



|Fp(m)| 2 , (12) 



with r„<,_y^ = 7.8 • 10" 6 MeV, T^ 77 = 4.6 ■ 10" 4 MeV 
and the form factors, 

F n o{p) = 1 + 6^0 p 2 , 



Lo =5. 5 GeV" 



(13) 



FM = 1 - 



A v = 0.676 GeV. (14) 



The above value of A v has been recently determined from 
the HADES data at 2.2 GeV beam energy [S3] and agrees 
reasonably well with the values found by NA60 [60] and 
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CB/TAPS [61]. As shown in [62], the dilepton decays of model of [68] • However, we note that recently a new form- 
the pseudoscalar mesons is expected to follow an anisotropic factor calculation has appeared [69] , whose results differ 
angular distribution, significantly from the ones given in |68j . 



dr, 



dcos6> 



oc l + cos 2 (6>), 



(15) 



where 9 is the angle of the electron momentum with re- 
spect to the dilepton momentum. This has been confirmed 
recently by HADES data [63] . All other decays are treated 
isotropically in our model. 
The parametrization of the w Dalitz decay, 



dr 



cj— ^Tr"e T e 



2a r, 



dfi 



3tt p 



M 



x li^oor, 



2 



^lf 



(t'i 



I^(m)I 2 



(a*-w + a*i* : 



3/2 

(16) 
(17) 



is adopted from [64 | [41 ] with r t 



UJ— >-7T u 7 



0.703 MeV, A a 



0.65 GeV and r u = 75 MeV. Here we note that the form 
factor of the u Dalitz decay is also well-constrained by 
data [50] . 
For the Z\-Dalitz decay, we use the parametrization from 



For the other baryonic resonances we don't explicitly in- 
clude a Dalitz decay, but evaluate their contributions to 
the dilepton spectrum through the two-step process R — > 
Np — > Ne + e~ . In the transport-typical manner we cut 
the corresponding diagrams, separating the production 
and decay vertices of the resonance and neglecting any 
phases and interferences. Below the 2tt threshold, the p 
meson width becomes very small because here only the 
electromagnetic decay width is active. This smallness of 
the width, however, is counteracted to some degree by 
the propagator of the vi rtua l photon that enhances small 
dilepton masses, see eq. ( 11 ). In an alternative treatment, 
in which the N* resonances undergo direct Dalitz de- 
cay, an electromagnetic form factor at the AW* 7* ver- 
tex would mimick the p propagator. These two methods 
are fully equivalent if the phase relations between the de- 
caying resonance and the dileptons can be neglected and 
a corresponding form factor is used (our model relies on 
the assumption of strict VMD). Any interaction of the 
p meson between its production and decay, leading to a 
broadening of the p spectral function, could be absorbed 
into a medium dependence of the form factor. 

Further we include pn-Bremsstrahlung in phase-space cor- 
rected soft-photon approximation |7UI72j . which can be 
written as 



d-^A-fJVe+e- 

d// 



2a 



r A ^ 



N~/* 



37^-^*' 
a (to^ + ttt-n) 2 
16 m^mff 



(18) 
[{m A +m N ) - p \ 



21 3/2 



x [{mA~m N y-p 2 Y"\F A {p)\ 



(19) 



where we neglect the electron mass. The electromagnetic 
N-Z\ transition form factor F A (p) is an issue of ongoing 
debate. Unlike the other semileptonic Dalitz decays, it is 
poorly constrained by data. At least at the real-photon 
point (p = 0) it is fixed by the decay width r A ^Nj sa 
0.66 MeV [37] to \F A (Q)\ = 3.03, and also in the space- 
like region this form factor is well-constrained by electron- 
scattering data on the nucleon. However, it is basically 
unknown in the time-like regime, which is being probed 
by the A Dalitz decay. 

Theoretical models for the N-A transition form factor usu- 
ally assume one or more VMD-inspired peaks in the time- 
like region [66, 67, 68,69J. However, the data in the space- 
like region does not provide sufficient constraints to fix the 
behavior in the time-like region. 

Moreover, a VMD-like A form factor would imply a cou- 
pling of the A to the p meson, which has never been ob- 
served directly and could only play a role far off the A 
pole, where its strength is completely unknown 1701 . 
In order to demonstrate the uncertainty connected to this 
form factor, we will in the following use as an example the 
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where M is the mass of the dilepton pair, q, E and fl are 
its momentum, energy and solid angle in the pn center-of- 
mass frame and s is the Mandelstam's variable. Further, 
777i is the mass of the charged particle (proton), 7772 is 
the mass of the neutral particle (neutron) and cr^™ is the 
elastic pn cross section. 

pp-Brcmsstrahlung can not be treated in this simple ap- 
proximation, since it involves a destructive interference 
between the graphs involved. Due to this interference it is 
much smaller than the pn-Bremsstrahlung and therefore 
is being neglected here. 

Further we note that also the Bremsstrahlung contribu- 
tion involves a form factor, i.e. the time-like nucleon form 
factor. Just as the A transition form factor, it is not 
well-constrained in the time-like region and is usually ne- 
glected, also in recent OBE models [73 ]I74|I75| . 
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5 Dilepton spectra from elementary N+N 
collisions 

After describing the basic ingredients of the model, we will 
now present simulated dilepton spectra for various elemen- 
tary reactions (p+p and d+p). The HADES collaboration 
has performed measurements of dilepton spectra from ele- 
mentary p+p reactions at the beam energies listed in table 

H 



p + p at 1 .25 GeV 



-Ekin 


V^ 


Plab 


Plcp 


max 
Plcp 


1.25 


2.42 


1.98 


0.05 


1.8 


2.20 


2.76 


2.99 


0.10 


2.0 


3.50 


3.17 


4.34 


0.08 


2.0 



Table 4. Kinematic conditions of the elementary collisions 
measured by HADES and corresponding cuts on the single 
lepton momenta (all numbers in GeV). 



For the beam energy of 1.25 AGeV, also dp collisions have 
been measured. All sinrulated spectra have been filtered 
with the HADES acceptance filter (HAFT, version 2.0) 
[751177] , in order to take care of the geometrical accep- 
tance and resolution of the detector. In addition, a dilep- 
ton opening angle cut of 8 ee > 9° is applied in all cases, 
as well as the single-lepton momentum cuts listed in table 
|4j matching the experimental analysis procedure. 



5.1 p + p at 1.25 GeV 

The lowest HADES energy, i? kin = 1.25 GeV, correspond- 
ing to y/s 1 sa 2.4 GeV, is just below the rj production 
threshold, and also for p mesons there is only a small 
sub-threshold contribution from the low-mass tail of the 
p spectral function. 

This means that the dilepton spectrum is dominated by 
the 7T° and A Dalitz decays. One should note that at this 
energy, almost all pions are produced via excitation and 
decay of the A resonance. 

Both of these Dalitz decays involve a transition form fac- 
tor. But while the form factor of the ir° Dalitz channel has 
been determined experimentally to a reasonable precision 
[58] . the electromagnetic transition form factor of the A 
Dalitz decay is basically unknown in the time-like region 
(cf. previous sect.). 

However, the dilepton spectrum at E\& n = 1.25 GeV is 
only mildly sensitive to this form factor, since the energy 
is not large enough to reach the VM pole-mass region. As 
fig. [2] shows, the simulation profits from including a form 
factor (shaded band) which exhibits a moderate rise in 
the time-like region of small q 2 , but it is not sensitive to 
the actual VMD peak of such a form factor. Here we have 
used the form factor from [55], but we have also verified 
that using a standard VMD form factor yields virtually 
the same results for this energy. 

It is interesting to note that other calculations achieve a 
good agreement with the HADES data for pp collisions 
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Fig. 2. (Color online) Dilepton mass spectrum for pp at 
1.25 GeV, in comparison to the data from |22] , The different 
contributions are indicated in the figure. The hatched area in- 
dicates the effect of the A form factor. 



at 1.25 GeV without including any form factor for the A 
Dalitz channel [75j (which might be partly due to the dif- 
ferent width parametrization used). 

Furthermore we note that the slight overshooting in the 
pion channel is apparently due to the higher resonances, 
which are produced only in phase-space approximation. 
The missing treatment of proper angular momentum dis- 
tributions seems to interfere with the HADES acceptance 
here, however it does not seem to be a problem at higher 
energies. We have verified that the discrepancy disappears 
if all pions are produced exclusively via A excitation, ne- 
glecting contributions from higher resonances. 



5.2 d + p at 1.25 GeV 



In addition to the proton beam, also a deuteron beam with 
a kinetic energy of 1.25 AGeV has been used by HADES. 
Here, a trigger on forward-going protons has been set up 
in order to select the (quasi-free) np collisions, which are 
only accessible in this way. 

Due to the motion of the bound nucleons in the deuteron, 
the energy of the NN collisions is smeared out here, com- 
pared to the proton-beam case, with a tail reaching above 
the ^-production threshold. The momentum distribution 
of the nucleons is determined by the deuteron potential, 
which in our simulations is given by the Argonne V18 po- 
tential [78] , 

Fig. [3] shows the dilepton invariant mass spectrum for this 
reaction. While the 7r°-Dalitz channel in the low-mass re- 
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d + p at 1 .25 GeV 




According to our analysis, the most probable candidate 
to fill the missing yield are indeed p-like contributions. 
The radiation from internal pion lines is one such graph 
which we miss; this channel mainly contributes at large 
masses [75 . Furthermore, the subthreshold p production 
via resonances could be underestimated on the neutron by 
our model. Analogous to the rj case, it might be enhanced 
over pp — > p°X. And finally, channels like np — > dp (re- 
lated to the radiative capture) could contribute, which are 
completely unknown. 

The discrepancy of data and theory for the d+p reac- 
tion is specific for this reaction at this particular energy; 
the results for nuclear collisions to be discussed later do 
not show such a disagreement. We note that the observed 
cross section represents only about 15 - 20% of the actual 
cross section; the rest is being cut away by the acceptance 
filter. Thus, any deficiencies, for example, in the angular 
distribution of our dileptons could show up in rather large 
errors of the spectra after the acceptance cuts have been 
performed. 



0.1 0.2 0.3 0.4 0.5 0.6 0.7 
dilepton mass m ee [GeV] 



Fig. 3. (Color online) Dilepton mass spectrum for d+p at 
1.25 GeV, in comparison to the data from |22| . 



gion shows a similarly good agreement as in the pp case, 
the data points at larger invariant masses are underesti- 
mated by a factor of two or more. 

A stronger A channel can apparently not explain the shoul- 
der in the data around 500 MeV, since it falls off too 
steeply, even when including a form factor. In addition 
to the enhanced r\ production in np — > rvprj, as described 
earlier, we have included a np — > dry channel, which dom- 
inates the rj production from np at threshold [45] . 
Unfortunately, the strong pn — >• pnp° channel is experi- 
mentally not so well known. In our model, the p° produc- 
tion in d+p at 1.25 GeV is dominated by the -Di3(1520) 
and S'ii(1535) resonances. The latter is enhanced in np 
(because of its dominant role in r\ production). The for- 
mer is assumed to be isospin-symmetric, which may not 
be the case. 

In an OBE-model study [75] it has been found that the 
radiation from internal pion lines (with the appropriate 
VMD form factor) gives a sizable contribution at large 
invariant masses. Such a diagram implicitly contains a p° 
propagator (through the form factor) , and gives additional 
p-like contributions on top of the resonance contributions 
included in our model. 

Moreover, we might underestimate the 'pure' Brcmsstrah- 
lung contributions, which do not involve resonance excita- 
tions, due to the soft-photon approximation. However, it 
is not expected that these terms would yield any dominant 
contributions [79,80,75 _. 

As recently argued in [81 j . the inclusion of a "radiative 
capture" channel np — >• de + e~, fixed via deuteron photo- 
disintegration, might give further contributions in the high- 
mass region. 



5.3 p + p at 3.5 GeV 

Fig. [4] shows a comparison plot of a GiBUU simulation to 
HADES data [24] for a proton beam of 3.5 GeV kinetic 
energy impinging on a fixed proton target. This is the 
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Fig. 4. (Color online) Dilepton mass spectrum for pp colli- 
sions at 3.5 GeV. Data from |24] , The hatched areas indicate 
the effects of the A form factor [68] and baryon-resonance con- 
tributions to the p production, respectively. The total is shown 
(from bottom to top) with A form factor (left- hatched), p res- 
onance contributions (right-hatched) and with both of these 
effects together. 
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Fig. 5. (Color online) Top: Resonance contributions to the 
p channel in the dilepton mass spectrum. Bottom: Resonance 
contributions to the p mass distribution. The dashed line in- 
dicates the vacuum pole mass of the p meson. For comparison 
we also show the p meson contribution from our earlier Pythia 
simulations 33 . 



highest beam energy (per nucleon) used by the HADES 
experiment and corresponds to a center-of-mass energy of 
y/s = 3.18 GeV. 

At this energy the r), uj and p production channels are fully 
open, and even the </> production becomes energetically 
possible. The data only show a hint of a <j) peak with very 



poor statistics, but it seems to be slightly underestimated 
by our simulation. 

The r\ and p production is dominated by the channels 
NN -» NNtttj and NN -> NNnp, respectively. In our 
model these are saturated by double-resonance excitation, 
cf. sec. [3] The 00 meson is presently produced in a non- 
resonant phase-space prescription through the exclusive 
and the ttoj channel. 

Under these assumptions, we get a very good agreement 
with the data over the whole mass range, as shown in 

fig. HI 

It is interesting to note that the shape of the p channel 
shown here differs significantly from the one obtained in 
our previous string- model investigations via Pythia |33) . 
which has been adopted for the PLUTO simulations in 
[24] . The latter is given by the lower dashed (green) line in 
Fig.|4l whereas the new resonance-model based treatment 
yields the upper dashed line. The p-shape effect is due to 
the production of p mesons via nucleon resonances, i.e. 
NN -> NB, -> NNp and NN -> AR -> NNnp, where 
the lighter resonances like e.g. Z?i3(1520) will preferen- 
tially contribute to the low-mass part of the p spectral 
function. Together with the 1/m 3 factor of the dilepton 
decay width, this results in a very flat distribution, which 
lacks a clear peak at the nominal mass, and dominates the 
dilepton spectrum in the intermediate mass region around 
500 - 700 MeV. 

The p spectral function is thus 'modified' already in the 
vacuum, simply due to the production mechanism via nu- 
cleon resonances. As seen in fig. [5] the p mass distribution 
in pp at 3.5 GeV peaks around 730 MeV, with an addi- 
tional shoulder around 500 MeV (due to low-mass reso- 
nances, mainly the Z?i3(1520)). This spectral shape is due 
to phase-space limitations and special resonance proper- 
ties. It differs significantly from the mass distribution re- 
sulting from a Pythia simulation [33], which lacks any 
resonance contributions. Similar effects were already ob- 
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Fig. 6. (Color online) Inclusive and exclusive A + production 
cross sections in different models (Fritiof 7.02, Pythia 6.4 
and the GiBUU resonance model), compared to data from |82j . 
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Fig. 7. (Color online) Transverse momentum and rapidity spectra of dilepton pairs from pp at 3.5 GeV in four mass bins. The 
hatched area indicates the effect of the A form factor. Data from I24I. 



served, e.g., in C+C reactions [53]. We stress here that this 
is not an 'in-medium' effect at all: It is solely caused by 
the production mechanism and occurs already in elemen- 
tary p+p collisions in the vacuum. This effect is crucial for 
understanding the intermediate mass region of the dilep- 
ton spectrum in pp collisions at 3.5 GeV (as seen in fig.[4| 
and might also play an important role at 2.2 GeV (see next 
section). 

The particular influence of the AT* (1520) resonance on 
dilepton spectra from NN collisions have already been in- 
vestigated in [51], where it was concluded that the N* (1520) 
can indeed give sizable contributions to the DLS and HA- 
DES spectra, but is subject to moderate uncertainties. 
It should be noted that the exact composition of the res- 
onance contributions to the p channel, and therefore also 
its exact shape, are not fixed by data so far, but rather 
represent an 'educated guess'. The resonance composition 
can be checked via ttN invariant mass spectra. 
Moreover, possible pA decay modes of certain resonances 
could give further contributions to the dilepton cocktail, 
as mentioned earlier. 

Comparing our cocktail to other transport models like 
HSD [55] or UrQMD [5J5] , one of the most significant dis- 
crepancies shows up in the size of the A channel. While 
in our model the A does not give any significant con- 
tribution to the total dilepton yield at E^in = 3.5 GeV 
(without a form factor), this is not so for the two other 
models. Both of them have a much stronger A channel, 
which even dominates the dilepton spectrum in the inter- 
mediate mass region around 600 MeV. We stress here that 



there are several factors of uncertainty in the A channel, 
for example the inclusive production cross section, but 
also the parametrization of the A decay width (hadronic 
as well as leptonic) and the completely unsettled question 
of the electromagnetic N-A transition form factor. 
Although the inclusive A production cross section is not 
that well known at £k; n = 3.5 GeV, one can get con- 
straints from the exclusive cross section, cf. fig. [6] as well 
as the inclusive one at lower energies (where it is fixed 
via pion production). Both constraints are respected in 
our resonance model, while e.g. the Fritiof model clearly 
overestimates the exclusive A + production, and in partic- 
ular does not seem to respect the correct isospin relations. 
On the question of the electromagnetic N-A transition 
form factor, it should be noted that in our simulations 
the Iachello model [68] agrees reasonably well with the 
data (depending on the contributions of other baryonic 
resonances), while a naive VMD form factor, as used e.g. in 
[87] . would clearly overshoot the data. 
In order to understand the underlying processes, it is not 
sufficient to consider only the mass spectrum. Other ob- 
servables can give further insight into the reaction dy- 
namics and can serve as a cross check for the validation of 
theoretical models. In order to compare to the data from 
[24] . we examine the transverse momentum and rapidity 
distributions in four different mass bins (see fig. [7]): 

• m < 150 MeV, dominated by the 7r° Dalitz channel, 

• 150 MeV < m < 470 MeV, dominated by the r] Dalitz 
decay, 
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Fig. 9. (Color online) Transverse momentum spectra of dilepton pairs from pp at 2.2 GeV in three mass bins. The hatched 
area indicates the effect of the A form factor. Data from I25I. 



• 470 MeV < m < 700 MeV, dominated by the direct p 
decay (possibly with contributions from the A Dalitz), 

• 700 MeV < m, dominated by the ui and p. 



Distinguishing several mass bins is useful in order to sep- 
arate the contributions of different channels. In all four 
mass bins, we achieve an excellent agreement with the 
HADES data 24J. In particular it should be noted that 
a stronger A channel would apparently destroy the very 
good agreement in the pt spectra, since it would yield 
too large high-pT contributions in the two mass bins of 
150-470 and 470-700 MeV. 



p + p at 2.2 GeV 
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Fig. 8. (Color online) Dilepton mass spectrum for pp at 
2.2 GeV. The hatched area indicates the effect of the A form 
factor. Data from 1251. 



5.4 p + p at 2.2 GeV 

A third, intermediate, beam energy of 2.2 GeV has been 
used for the HADES experiment. This energy is well above 
the rj production threshold and is just high enough to 
reach the pole mass of the light vector mesons, p and u, 
which dominate the high-mass part of the dilepton spec- 
trum (as seen in fig. fSp . The A channel plays a less im- 
portant role here, since it is buried underneath the strong 
•q and p channels. The u only gives a small contribution, 
since the energy is only just at the threshold of to produc- 
tion. 
The p channel exhibits slightly more structure here than at 

3.5 GeV, showing a moderate step around 550 MeV. This 
step marks the border between a low-mass part, which is 
dominated by the Di3(1520) resonance, and a high-mass 
part dominated by the Pi3(1720). In fig. [8] we show the 
contributions of these two resonances to the p channel, but 
omit the subdominant contributions of other resonances 
(for the sake of readability) . The resonance contributions 
indeed improve the agreement with the data, compared 
to the PLUTO cocktail, which only includes phase-space 
population of the p |25j . However, there are still minor 
deviations, which seem to suggest an underestimation of 
the Di3(1520) and an overestimation of the Pi3(1720) in 
our resonance cocktail at this energy. 

In fig. [9] we show the pt spectra for three different mass 
bins in comparison to the data from [25] . Our simulations 
give a better agreement with the data than the PLUTO 
cocktail shown in [35] in all three mass bins. Most notably, 
we get an improvement from the larger p contribution in 
the highest mass bin. 



5.5 Comparison to elementary DLS data 

In addition to the recently measured HADES data, also 
the elementary data measured previously by the DLS col- 
laboration are available for comparison with our model 
[88] , Unfortunately they are of inferior quality in terms 
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DLS: p+d 
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Fig. 10. (Color online) Dilepton mass spectra in comparison to DLS data [55]. Top: p+p, bottom: p+d. 
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of statistics and acceptance. However, more beam ener- 
gies have been measured than in the case of HADES, so 
that they can still provide additional contraints, which are 
useful for understanding the elementary cocktail. 
In order to compare to the DLS data, the GiBUU dilepton 
events have been filtered throug the DLS acceptance fil- 
ter, version 4.1, as available from [89| . In addition to the 
acceptance filtering, the events have been smeared with 
a Gaussian of width a — 0.1m co , in order to account for 
the mass resolution of the detector. No further cuts have 
been applied. The kinematics of the reactions measured 
by DLS are summarized in table [5] At each of the given 
energies, a p+p and p+d reaction was measured. 



p + Nb at 3.5 GeV 
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4.88 


3.56 


5.74 



Table 5. Kinematic conditions of the elementary collisions 
measured by DLS (in GeV). 



The comparison of the GiBUU model results to the DLS 



data is shown in Fig. 10 As before, we show the effect of 
the A transition form factor as a hatched band and note 
that it slightly improves the agreement with the data in 
almost all cases. 

Is is apparent that at the medium beam energies there 
is a reasonable agreement, both in p+p and p+d. The 
largest deviations are visible at the highest beam energy 
of 4.88 GeV, which is already at the border of validity of 
our resonance model. Apparently the inclusive production 
of p and w mesons is underestimated there. 
The underestimation at the lowest energy of 1.04 GeV is 
similar to that seen in the HADES experiment at a com- 
parable energy (see Fig. [3]). Since at this somewhat lower 
energy the r\ production plays no role, the discrepancy 
seems to indicate a problem with the A or Bremsstrah- 
lung contributions. However, we note again that the popu- 
lation of the A resonance is constrained rather well by the 
pion and total cross sections, which we describe rather 
well (see Fig. fl]). On the other hand, the decay of the 
resonance is fixed by the electromagnetic coupling at the 
photon point, so that there is no ambiguity there. Further, 
form factors have only little influence at such low eneries. 
We thus have to conclude that we have no explanation for 
the discrepancy yet and note that related, earlier calcula- 
tions similarly underestimated the DLS dilepton yield at 
this lowest energy [19]. 
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Fig. 11. (Color online) Dilepton mass spectrum for p+Nb at 
3.5 GeV, showing all contributing channels with vacuum spec- 
tral functions. The hatched area indicates the effect of the A 
form factor. Preliminary data taken from [90], scaled to fit the 
7T° and r\ yield. 



GeV, we filter our dilepton events through the HADES 
acceptance filter and cut on 0.08 GeV < p\ cp < 2.0 GeV 
and 8 CC > 9° . The level of agreement is similar to the p+p 
reaction at the same energy. Note, however, that the data 
are not absolutely normalized in terms of a cross section 
yet. Therefore we have scaled the data points to match the 
simulation in the low-mass region, which is dominated by 
the 7T° and r\ Dalitz channels. Moreover, the data have not 
been fully corrected for all detector effects yet, which is 
the reason for a slight shift of the mass scale (on the order 
of 1%), which is visible at the w peak [ST]. It is evident 
that the data can be quite well described if the electro- 
magnetic A decay width does not contain the form factor 
of ref. [68] which would create a hump in the spectrum 
around 0.6 GeV. 

In contrast to the NN collisions in the preceding chapter, 
here we neglect the p-meson contributions below the 2ir 



threshold (due to numerical reasons). As seen in fig. 11 



they do not contribute significantly to the total dilepton 
yield. 



6.1 In-medium effects 



6 Dilepton spectra from p+Nb collisions 

Fig. [TT] shows simulated dilepton spectra for p+Nb colli- 
sions at 3.5 GeV using vacuum spectral functions, com- 
pared to preliminary data from [90] . As for p+p at 3.5 



In p+Nb reactions there are additional effects, compared 
to the elementary p+p reactions. First of all, the primary 
p+N collisions will be nearly identical, apart from bin- 
ding effects and some Fermi smearing and Pauli blocking, 
but besides p+p also p+n collisions play a role. Further- 
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Fig. 12. (Color online) Dilepton mass spectra for p+Nb at 3.5 GeV. Comparison of different in-medium scenarios (vacuum 
spectral functions for the vector mesons, collisional broadening, 16% mass shift, collisional broadening plus mass shift). Left: 
Total spectrum, center: p contribution, right: w contribution. Preliminary data taken from [90], scaled to fit the ty° and r\ yield. 



more, the produced particles undergo final-state interac- 
tions within the Nb nucleus, and processes like meson ab- 
sorption and regeneration may become important. The 
secondary collisions will on average have lower energies 
than the primary N+N collisions. Finally also the vector- 
meson spectral functions may be modified in the nuclear 
medium. 

The propagation of particles with density-dependent spec- 
tral functions (usually referred to as "off-shell propaga- 
tion") poses a particular challenge. Our approach to this 
problem is based on the off-shell equations of motion of 
test particles, as given in [92] and [93]. Such an off-shell 
treatment is necessary for including in-medium modifica- 
tions of the spectral functions (e.g. collisional broadening 
of the vector mesons). The collisional width inside a nu- 
clear medium of density, p, can be related to the collision 
cross section, onx, in low-density approximation as 



r. 



coll = P {VrdCTNX) 



(24) 

where v TC \ is the relative velocity and the brackets indi- 
cate an integration over the Fermi momentum of the nu- 
cleons. This collisional width will in general depend on the 
momentum of the involved particle, X. In order to avoid 
numerical difficulties connected with the appearance of 
superluminous test particles, we neglect the momentum 
dependence and use the simplified form, 



r„ 



r p 
Po 



(25) 



where po — 0.168 fm~ is the normal nuclear matter den- 
sity. The value of To should on average match the momen- 
tum-dependent width as obtained from the collision term. 
We typically use P = 150 MeV for the p and r = 80 MeV 
for the to meson. More details on off-shell propagation in 
the GiBUU model in general can be found in [26] . 
The mass spectrum above 500 MeV can receive modi- 
fications from the inclusion of in-medium effects in the 



in-medium scenarios: The first one includes a collision- 
ally broadened in-medium width, while the second one 
assumes a pole-mass shift according to 



m*(p) = mo I 1 — a 



P_ 
Po 



(26) 



vector-meson spectral functions. Fig. 12 shows the typical 



with a scaling parameter, a = 16%. The third scenario 
combines both of these effects. The modifications intro- 
duced by these scenarios are roughly on the same order of 
magnitude as the systematic errors of the data, and so far 
there is no clear evidence for medium modifications of the 
vector-meson properties in cold nuclear matter from the 
HADES data. However, it looks as if a mass shift tends to 
deteriorate the agreement with the data. 
Regarding the w absorption, it should be noted that the 
GiBUU implementation yields an average collisional width 
of roughly r = 80 MeV. This appeared too low to explain 
the transparency-ratio measurement of [94] , which seemed 
to demand values of 130 to 150 MeV. For the HADES 
dilepton data, such a discrepancy currently does not seem 
to exist. 

However, one should keep in mind that a statement about 
to absorption depends on a number of prerequisites. For 
example, one needs to have the p contribution well under 
control, since it represents a large background under the u> 
peak. Given the discussion about resonance contributions 
to the elementary p production, this is already not a trivial 
task, even more complicated by possible in-medium mod- 
ifications of the p meson. Furthermore, the size of the uj 
peak in pNb crucially depends not only on the production 
cross section in pp collisions (which is well determined via 
the elementary pp data at 3.5 GeV), but also in pn, which 
is unknown. We assume u production cross sections which 
are isospin-independent, i.e. equal in pp and pn. 
In addition to the in-medium modifications of the vec- 
tor mesons, also the baryonic resonances can receive sim- 
ilar modifications in the medium. Since the production 
via baryon resonances is particularly important for the p 
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meson, in-medium modifications of these resonances can 

lead to further modifications of the p contribution to the 

dilepton spectrum, which should be considered in future 

investigations. 

The pt and rapidity spectra for p+Nb are depicted in 

fig. 13 with the same mass binning as in the p+p case. 



The shown px and rapidity spectra do not include any 
in-medium effects for the vector mesons and are not sig- 
nificantly sensitive to such modifications. 
For a further discussion of the nuclear effects, it is useful 
to consider the quantity 



After fixing the model with the constraints given by the 
elementary N+N collisions, the p+Nb reaction at 3.5 GeV 
is reasonably well described by the GiBUU transport model, 
using the same input and without requring any in-medium 
mass shifts. According to our model, the p+Nb data show 
only a limited sensitivity to collisional broadening of the 
p meson. 

These results also provide the basis for a further investi- 
gation of the heavy-ion collisions at SIS energies measured 
by the HADES collaboration pJO p T p g ] ■ 



-RpNb — 



r pNb-j-c+o-X 
Opp^c+c-X 



°"pp->x 

CpNb^X 



(27) Acknowledgments 



i.e., the ratio of dilepton yields in pNb vs. pp, normalized 
to the total cross section for these reactions (whose ratio 
is roughly cr p Nb-s.x/c P p^x ~ 25.0 in our simulations). If 
medium effects are negligible, this quantity will be unity. 
Therefore, any deviation from unity indicates medium ef- 
fects such as, e.g., absorption (R < 1) or secondary pro- 
duction (R > 1) 



Fig. 
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shows -RpNb as a function of the 
dilepton momentum in four different invariant-mass bins, 
with the contributions from the different source channels. 
While -RpNb is relatively flat in the ir° region, the higher 
mass bins show a strong enhancement at low momenta, 
which can be understood as secondary particle production 
and/or elastic rescattering. The high momentum region in 
all mass bins tends to show a slight depletion, connected 
to absorption. 

The observable i? P Nb could also help to pin down the rel- 
ative contribution of the A Dalitz channel to the dilepton 
spectrum. As can be seen in Fig.[l4j the A's ratio is rather 
large, due to the enhanced production of the A + '° charge 
states in pn collisions, relative to pp. The isospin factors 
for NN — > NA + '° are a factor of two larger in pn than in 
pp. 

This isospin dependence could provide additional constraints 
for distinguishing the p and A contributions in the inter- 
mediate mass range of 470 - 700 MeV. Since the p channel 
dominates our simulated cocktail in this mass range (with- 
out a A form factor) , the total value of i? P Nb roughly fol- 
lows the i?-value of the p channel. If the spectrum would 
be dominated by the A Dalitz channel in this mass range, 
then the total value of i? P Nb would be more similar to the 
A's R- value. 



7 Conclusions 

We have shown that the HADES data from elementary 
N+N collisions can be described consistently by an ex- 
tended resonance model over the whole range of beam en- 
ergies. We have set up such a model based on the earlier 
resonance model approach by Teis et al. 
For describing the dilepton mass spectrum at the high- 
est beam energy of 3.5 GeV, an essential ingredient is a p 
spectral function, which is modified through the produc- 
tion via nucleon resonances, with an enhanced low-mass 
contribution from low-lying resonances like the £>i3(1520). 
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